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Abstract 

We describe a few features of the Stern sequence and of a closely 
related sequence obtained by adding a sign-twist in the recursive defi- 
nition of the Stern sequence o 

1 Main results 



o. 

v-^ ' In this paper, we identify a complex sequence (a n ) n6 N with the correspond- 

SZ ■ ing function a : N i — > C We write thus always a(n) instead of a n . 

The Stern sequence or Stern-Brocot sequence with first terms given by 

0,1,1,2,1,3,2,3,1,4,3,5,2,5,3,4,1,5,4,7,3,8,5,7,2,7,5,8,3,... 

(cf. sequence A2487 of [9j) is the integral sequence s : N — > N recursively 
F*"- i defined by s(0) = 0, s(l) = 1 and s(2n) = s(n),s(2n + 1) = s(n) + s(n + 1) 

for n > 1. It is closely related to the Farey tree and induces a one-to-one 
ly-) ■ map n i — > s(n)/s(n + 1) between N and non-negative rational numbers, cf. 

[5] or Chapter 16 of [TJ. It is also an example of a 2— regular sequence, see 
Chapter 16 of [2]. 

The following result gives a different, perhaps not very well-known, de- 
scription of the Stern sequence. 

S^ Proposition 1.1 s(n) equals the number of distinct subsequences of the 

form, 1, 101, 10101, . . . = {1(01)*} in the binary expansion e\ . . . e\e§ of n = 
Ei=o^2 fc (where eo,...,^^ {0,1}). 

Froposition 11.11 is in fact a particular case of Froposition 12.31 an easy 
result concerning rational series in non-commuting variables. 

Example The binary expansion 1011 of 11 = 2 3 + 2 1 + 2° contains the 
following five subsequences (highlighted by bold letters) 

1011 , 1011 , 1011,1011 ,1011 



1 Keywords: Stern sequence, automatic sequence, regular sequence. Math, class: 11B85 



of the form 1(01)*. We have a(ll) = s(5) + s(6) = s(2) + s(3) + s(3) = 
«(l)+2(a(l) + a(2)) = 5«(l) = 5. 

Proposition 11.11 allows to parametrize Stern sequences by counting a 
subsequence of the form l(01) fc with weight w k , see Proposition 12.11 for 
formulae. For n = 11 we get for instance 3 + 2w. 

In this paper we introduce a related sequence which will be called the 
twisted Stern sequence since it is obtained by twisting the recursive definition 
of the Stern sequence with a sign. More precisely, we define the twisted 
Stern sequence £(0), t(l), . . . recursively by t(0) = 0, t(l) = 1 and t(2n) = 
—t{n),t{2n + 1) = — t(n) — t(n + 1) for n > 1. It starts as 

0,1,-1,0,1,1,0,-1,-1,-2,-1,-1,0,1,1,2,1,3,2,3,1,2,1,1,0,-1,... . 

An inspection of these first few terms shows already some striking similarities 
between the Stern sequence and its twisted relative. The aim of this paper 
is to describe a few properties of the Stern sequence and its twist. 

The following result (the identity for s(n) in assertion (i) is probably well- 
known to the experts) is an illustration of the similarities between these two 
sequences: 

Theorem 1.2 (i) We have 

s(2 e +n) = s(2 e -n) + s(n) 

t(2 e + n) = (-l) e (s(2 e - n) - s{n)) 

for all e > and for all n such that < n < 2 e . 
(ii) We have 

t(3 -2 e + n)= t(6 • 2 e - n) = {-lfs{n) 

for all e > and for all n such that < n < 2 e+1 . 

The failure for n > 2 e+1 of the formula 

t(3-2 e + n) = (-l) e s(n) 

given by assertion (ii) can perhaps be mended by the following conjectural 
identity based on experimental observations. 

Conjecture 1.3 There exists an integral sequence u(0),u(l),u(2), . . . such 
that we have 



J2 «(3 • 2 e + n)z n = (-l) e J2 u ( n > n ' 2e E 

n=0 \n=0 / \m=0 

for all e G N. 



s(m)z m 



If the conjecture holds, the ordinary generating function of the sequence 
u(0),u(l), ... is given by 



oo 



X>( 



Y,n=O s ( n ) zn 



n=0 

and it starts as 

1 - 2z 2 - 2z 5 + Az G + 2z 7 - 6z s + 4z 9 + 2z 10 - Qz 11 + 8z 12 + ... 

The first equality in assertion (ii) of Theorem 11.21 shows that the finite 
sequences (-l) e t(3 • 2 e ), (-l) e t(3 ■ 2 e + 1, . . . , (-l) e t(6 ■ 2 e ) of length 3 e + 1 
are palindromic sequences of natural integers. The first few such sequences 
are 

110 

112 110 

0112132312110 
0112132314352534132312110 

with boldfaced l's at one third and two thirds highlighting the underlying 
partial self-similarity structure. The subsequence lying between the two 
boldfaced l's appears also at the beginning of [8j. One notices that all 
sequences start and end with zero and that all existing central elements are 
equal to 2. 

The polynomials defined by these palindromic sequences are described 
by the following result: 

Theorem 1.4 The polynomials 

3-2 e 



^ e = (-l) e Y J t(3-2 e + n)z r ' 



n=0 

have the factorisations 

e-2 



i> e = Z (l + Z ^)(l + Z + Z lrH(l- Z l" +Z l » +1 )e-l-n (1) 

n=0 

= z(i+z 2e )n(i+z 2 ' i +z 2 " +i ). (2) 



n=0 

Remark 1.5 Theorem \l-4\ implies the identity 

oo od e—1 



Y,t(n)z n = z-z 2 + ^(-l) e z 3 - 2e+1 (l + z 2e ) IJ(1 + z 2 " + z 2n+1 ) . 

n=0 e=0 ra=0 



Assertion (ii) of Theorem 11.21 yields ]hn e - +00 ip e = J2^=o s ( n )t n - The 
factorisation ([2]) of Theorem 11.41 gives a new proof of the following result 
due to Carlitz (see [6]): 



Corollary 1.6 We have £~ =0 s(n)z n = zU™ =Q {l + ^ + z2 " +1 )- 

(A direct proof of Corollary 11.61 is straightforward: The series U(z) de- 
fined by the right-hand-side starts as z + . . . = s(0) + s(l)z + . . . and its even, 
respectively odd, subseries are given by U(z 2 ), respectively Q + z) U(z 2 ). 
Its coefficients satisfy thus the same recursion relations as the elements of 
the Stern sequence.) 

The Carlitz factorisation of Corollary 11.61 implies that s ^2m = Qs{n)z n has 
no non-zero roots in the open unit disc. This is not true for the ordinary 
generating series Y^=o^{ n ) zn of the twisted Stern sequence which has (in- 
finitely?) many non-zero roots in the open unit disc. 

Given a natural integer k > 2 and a natural integer i, we consider the 
endomorphism p(i) of the vector-space (or module) of formal power series 
defined by 

(oo \ oo 

J2a(n)z n \ = ^a{i + nk)z n . 
n=0 / n=0 

The k— kernel of a formal power series A is the smallest vector space (or mod- 
ule when working over a ring) V containing A such that p(0)V, . . . , p{k — 
1)V C V. A formal power series is k— regular if its /c— kernel is finitely gen- 
erated. Easy examples of k— regular series are polynomials and ordinary 
generating series of periodic sequences, k— regular power series form a vec- 
tor space (or module) which is preserved by many natural operations such 
as derivation, product, Hadamard product, shuffle product, ..., see [2] for 
details. The set of k— regular sequences with coefficients contained in a fi- 
nite set (eg. in a finite field) coincides with the set of so-called A;— automatic 
sequences, see Theorem 16.1.5 of [2J. 

A sequence a(0), a(l), . . . is called /c— regular if its ordinary generating 
series Y^=o a i n ) zn ls a &— regular formal power series. 

The following result is also a consequence of the Carlitz factorisation: 

Theorem 1.7 The logarithmic derivative 

H(z) = ± log (± s(n + l)zA = ^^ ^T 

of YlnLo s ( n + I)- 2 " * s 2— regular. More precisely, H(z) is defined by the 
functional equation 



Coefficients of H(z) appear as sequence A163659 in [9|. 
2— regularity of the logarithmic derivation H{z) is a special case of the 
following result, perhaps already known to Schiitzenberger: 

Theorem 1.8 Given d k— regular series A\(z), . . . ,A d {z) over some com- 
mutative ring R, d linear forms L\{x\, . . . , xj), ■ ■ ■ , L d (xi, . . . , x d ) in d un- 
knowns with coefficients in R[z] and d constants a±, . . . ,a d € R such that 
a% = Ai(0) + Li(ai, ■ ■ ■ , a d ) (mod z) for i = 1, . . . , d, the system of equa- 
tions 

U 1 (z) = A 1 (z)+L 1 (U 1 (z k ),...,U d (z k )), 

U d {z) = A n (z) + L n {U x {z k ), ..., U d (z k )) 

determines a unique set of d k— regular sequences Ui(z), . . . , U d (z) with con- 
stant coefficients a.% = Ui(0) for i = 1, . . . , d. 

Remark 1.9 We have ol% = Ai(0) if the linear form Li has all its coeffi- 
cients in zR[z]. 



The series 



oo _. 

a ^ = n r^^ 



n=0 

satisfying A{z) = {Y^=q z n ) A(z 2 ) is not 2— regular (see Remark 11.111 be- 
low). This shows that Theorem 11.81 can not be extended to equations with 
linear forms having k— regular series as coefficients. 
Theorem 11.71 can be generalised as follows: 

Theorem 1.10 Let P(z) be a polynomial with constant coefficient 1. Then 
the series 

oo 

A=\\P{z kn ) 

ra=0 

is k— regular. 

Moreover, if all roots of P{z) are complex roots of 1 having finite order, 
then the logarithmic derivative B = A' /A of A is also k— regular. 

Remark 1.11 Given a k— regular series A(z) with constant coefficient 1, 
the product n^Lo A(z kr ' ') is generally not k— regular. Indeed, starting with 
the 2— regular series A(z) = l + z + z 2 + ... = j^, the coefficient b(n) in 
the series 



B = Y J b(n)z n = f[ T ^ 

n=0 n=0 



OH 

Z A 



counts the number of partitions of n into powers of 2, see sequence A123 in 
J2I, and log(b(2n)) is asymptotically equal to 2 ^ 2 v°9 Toqn) ( see eaua ^ on 
1.3 in %4%) which is incompatible with k— regularity of B by Theorem 16.3.1 
in Jjf . 



12 13 2 


3 


1 


4 


3 


-10 110 


-1 


-1 


-2 


-1 



Remark 1.12 Another famous sequence illustrating Theorem \1.10\ is the 

sequence 

oo oo 

l[(l-z 2n ) = Y / (-l) tm ^z n 

n=0 n=0 

related to the Thue-Morse sequence n i — > tm(n) defined by digit-sums mod- 
ulo 2 for binary expansions of natural integers. 

Another link between the two sequences s and t is given by determinants 
of 2 x 2— matrices. For n > 1 we consider the matrix 

M (n)=( S{n) S{U + 1) 
M[n) \t(n) t(n + l) 

with first row two consecutive terms of s and and second row the two corre- 
sponding consecutive terms of t. The first matrices are 2 x 2— submatrices 
defined by two consecutive rows of 

525351547 
-10 112 13 2 3 

Injectivity of the map n i — > s(n + l)/s(n) shows that the matrices 
M{n) are all different. Their determinants are characterised by the following 
result. 

Theorem 1.13 We have | det(M(n))| = 2 for all n > 1. More precisely, 

det(M(n)) = -2(-l) fc if 2 k < n < 2 k+1 . 

In particular, all matrices M{n) are invertible for n > 1. 

Since the recursive definitions of the integral sequences s and t differ 
only by signs (and since they satisfy the same initial conditions) they have 
the same reduction modulo 2 characterised by the following (easy) result, 
already contained in |10j : 

Proposition 1.14 The integers s{n) and t(n) are even if and only if n is 
divisible by 3. 

The determinants of the non-singular matrices M(n) are thus in some 
sense as small as possible: Indeed, since the sequences s and t coincide 
modulo 2, a matrix M(n) involves either a column consisting of even in- 
tegers or all its four entries are odd integers and such matrices have even 
determinants. 

We end this paper with a last result, going back to Stern (see [IU]) for 
the first part of assertion (i): 



Theorem 1.15 (i) The integer s(n) divides s(n — 1) + s(n + 1) for every 
n > 1. More precisely, we have 

s(n — 1) + sin + 1) 

-i VV = 1 + 2«2(n 

where the 2— valuation V2(n) is defined as the exponent of the highest power 
of 2 dividing n. 
The function 



c(z) = y s(n - 1) + g(ra+1) 

satisfies C(0) = and 



s(n) 

71=1 



C{z) = z\±^ + C(z') 

and is 2— regular. 

(ii) The integers t(n) and t(n — 1) + t(n + 1) are both zero if n = 3 • 2 k . 
They are both non-zero otherwise and t(n) divides t{n — 1) + tin + 1). More 
precisely, we have 

^- l ^l {n+l) =l + 2v 2 (n) 
t{n) 

if n (£ {2 n , 3 • 2 N } ; (t(0) + t(2))/t(l) = -1 and 

t(2 e - 1) + t(2 e + 1) 



t(2«) 

/or aZZ e sitc/i i/iai e > 1. 



1 + 2(e - 2) 



The rest of this paper is organized as follows: The next section contains 
a proof of Proposition 11.11 and a few complements. 

Section [3] is devoted to the (easy) proofs of Theorem 11.21 and to a few 
more formulae and conjectures involving the Stern sequence and its twist. 

Section[J]contains the proofs of Theorems ll.4^1.7U1.8iri.l0l and Corollary 
Ol 

Section [5] contains the easy proof of Theorem 11.131 and a few related 
results. 

Section [6] consists of the short proof of Theorem 11.151 

2 Proof of Proposition 11.11 and a few comments 

We give first a bijective proof of Proposition 11.11 We describe then briefly 
a weighted version of the Stern sequence counting subsequences of the form 



1(01)* with weights encoding their length. We give also a generalisation of 
Proposition 11.11 therefore providing a (sketch of a) second proof for Propo- 
sition 11.11 

Proof of Proposition [TTT1 We call a subsequence of a binary expansion 
B{n) admissible if it is of the form 1(01)*. For example, -B(ll) = 1011 = 
e3£2eieo has five admissible subsequences given by the set 

{e 3 ,e 1 ,eo,e3e2ei,e3e 2 eo} . 

Since the number b(n) of such subsequences satisfies clearly b(2n) = b(n), 
the equality s{2n) = s(n) shows that we can restrict our attention to n odd. 
We consider the two cases n = 4ra + 1 and n = An — 1. 

If w is an admissible subsequence of B(An + 1), then the digit t\ = 
of the binary expansion of 4re + 1 is either contained in w or not. In the 
first case, admissibility of w shows that w contains also the last digit eo 
of B(4n + 1) and removal of eieo = 01 from w yields a bijection between 
such admissible subsequences and admissible subsequences of B(2n). In 
the second case where e\ is not involved in w, we get a bijection between 
such admissible subsequences and admissible subsequences of the binary 
expansion B{2n + 1) = ... e3£2£o of the integer 2n + 1. This shows the 
identity b(4n + 1) = b(2n) + b(2n + 1). 

We consider finally the case of an admissible subsequence w of B(4n—1). 
If eo is not contained in w, then w can be associated with an admissible sub- 
sequence of B{4n — 2) or equivalenty of B(2n — 1). Denoting by I the least in- 
teger such that B(4n—1) = aOl , we consider now an admissible subsequence 
w of B(4n — 1) which contains eo- If the admissible subsequence w is not of 
the form /3e;eo, we transform it into the admissible subsequence /3ii = /31 of 
B(4n) = aiiei-i . . . eo = alO' or equivalently of B(2n) = alO i_1 (obtained 
from B(4n) by erasing the last digit never involved in an admissible sub- 
sequence). If w = /3e/eo we transform it into the admissible subsequence (3 
of B{4n) or equivalently of B(2n). This shows b{4n — 1) = b{2n — 1) + b(2n) 
and ends the proof. □ 

2.1 A weighted variation of the Stern sequence 

We denote by S(n) G N[w] the weighted number of subsequences of the 
form 1(01)* in the binary expansion B(n) of n, giving the weigth w k to a 
subsequence of the form l(01) fe . Similarly, we introduce S e (n) G N[w] as the 
weighted number of subsequences of the form (10)* in the binary expansion 
B(n) of n, with weigth w for a subsequence of the form (10) . 

Proposition 2.1 (i) Evaluating the polynomial S{n) G N[w] atw = 1 yields 
the Stern sequence. 



(ii) The sequences S(n) and S e (n) are uniquely determined by the initial 
conditions 5(0) = 0,5(1) = 5 e (0) = 5 e (l) = 1 and the recursive formulae 



5(2n) 


= S(n) 


5(2n + l) 


= S(n) + S e (n) 


S e (2n) 


= wS(n) + S e (n 


5 e (2n + 1) 


= 5 e (n) . 



(Hi) The sequence S(n) is also uniquely determined by the initial con- 
ditions 5(0) = 0, 5(1) = 1 and by the recursive formulae 

5(2n) = S{n) 
5(4n + l) = wS(2n) + S{2n + l) 
5(4n-l) = 5(2n - 1) + 5(2m + 1) + (w - l)5(2m) 

where n = 2 a (2m + 1). 

Remark 2.2 Klavzar, Milutinovic and Petri have studied a different family 
of polynomials closely related to the Stern sequence by considering Bq = 
0, B\ = 1, Bin = tB(n) and i?2n+i = B n + B n+ \, see Q/ for details. 

Proof of Proposition 12.11 Assertion (i) is obvious. 

In the sequel, we use the notation introduced above during the proof of 
Proposition 11.11 

The initial values for S(n) and S e (n) in assertion (ii) are easy to check. 
The identity 5(2n) = S(n) is obvious since admissible subsequences of B(2ri) 
never involve the last digit eo = in the binary expansion B(2n) of 2n. 

Admissible subsequences of B(2n + 1) not containing the last digit eo of 
B{2n + 1) are in weight-preserving bijection with admissible subsequences 
of B{2n) or of B(n). Removal of eo induces a weight-preserving bijection 
between admissible subsequences of B(2n + 1) involving the last digit eo of 
B{2n + 1) and monomial contributions to S e {n). This proves 5(2n + 1) = 
S(n) + S e (n). 

Monomial contributions to 5 e (2n) not involving the last digit eo of B(2n) 
are in (weight-preserving) bijection with monomial contributions to S e (n). 
Removing the last digit of monomial contributions to 5 e (2n) involving the 
last digit eo of B{2n) yields admissible subsequences of B(n) with weight 
reduced by 1. This shows 5 e (2n) = S e {n) + wS e {n). 

The identity S e (2n + 1) = S e (n) is due to the fact that monomial con- 
tributions to 5 e (2n + 1) never involve the last digit eo = 1 of B(2n + 1). 

Assertion (iii) follows from the bijections used in the proof of Proposition 
11.11 We leave the details to the reader. □ 



2.2 Counting weighted subsequences and subfactors 

A famous result by Schiitzenberger implies essentially an identification of 
k— regular sequences with the set of rational formal power series in k non- 
commuting variables. (One has to be a little careful with leading zeros. A 
way of dealing with them is to consider only formal power series involving 
no monomials starting with the variable xq associated to the digit 0.) 

Proposition 11.11 is then a particular case of the following well-known 
result which we give without proof. (A proof of stronger statements can be 
found in [3].) 

Proposition 2.3 Let Abe a rational formal power series in k non-commuting 
variables. Then the shuffle product of A with 1 _, * * and the ordinary 
non- commutative product -. — -, — \ — ; — ^A-, — -, — \ — ; — r are both rational. 

The shuffle-product counts subsequences encoded and weighted by A in 
k— ary expansions of natural integers and the ordinary product in Proposi- 
tion [2]3] counts subfactors (encoded by A) in A;— ary expansions. 

Proposition 11.11 corresponds to the case where 

1 

A = X\ = X\ + X\XqX\ + X\XqX\XqX\ + X\X§X\XqX\ + . . . 

1 — XqX\ 

respectively 

. W 2 

A = X\ = X\ + WX\XqX\ + W XlXo^lXoXl + . . . 

1 — XqXi 

in the weighted case with Xq,x\ non-commuting variables and w a central 
variable. 

A famous example counting subsequences or subfactors (reduced to 1) 
is given by the Thue- Morse sequence corresponding to A = x\. 

Another famous example counting subfactors is given by the Rudin- 
Shapiro sequence associated to A = x\ . 

3 Proof of Theorem 11.21 and more formulae 

3.1 Proof of Theorem IT21 

Proof of assertion (i) For n = we have 

a(l + 0) = s(l) + s(0) = 1 + = 1, a(l + 1) = s(0) + s(l) = + 1 = 1 

and 

t(l + 0) = s(l) - s(0) = 1-0=1, t(l + 1) = s(0) - s(l) = - 1 = -1 . 



10 



The proof is now by induction on e. If n is even we have 

S (2 e + n) = s{2^ + 2) = S (2 e -! - 2) + S (|) = s(2 * - „) + s(n) 



and 

t(2 e + n) = -t(2^ 1 + 2) 

= -(-ir 1 ^ 6 - 1 -^-^)) 

= (-l)^( s (2 e -n)- S (n)) 
If n is odd, we have 

s(2 e + n) = s(2 e " 1 + 2-^1) + s(2 e " 1 + ^±i) 

, . i n — 1. .n — 1. . . i n + l N ,n + 1. 

= s (2 e - 1 -^-) + s (^-) + s (2 e - 1 -^-) + S (^-) 

= s(2 e -n) + s(n) 
and 

t(2 e + n) 
= -t(2 e " 1 + )-t(2 e ~ 1 + — — ) 

= (-ir (.(2-' - S^i) - .(^) + .(2- - »±i) - .(2±i) 

= (-l) e W2 e -«)-»(»)) 

Proof of assertion (ii) These formulae are easy to establish for e = 0. 
For even n we have 

t (3 • 2 e + n) = -t(3 • 2- 1 + 2) = _(_i)e-i s( ^ = (-l)^(n) 

and 

t(6 . 2 e - „) = -t(6 • 2 6 " 1 - 2) = _(_i )e -i s( |) = (-l) e ,(n) 

and for odd n we get 

t(3-2 e + n) = -t(3 • 2 6 " 1 + ^-ti) - t(3 • 2 6 " 1 + ^-=^) 

- _(-i)-i( fl (!L±i) + a (!izJi)) = (-i)« 8 ( n ) 



and 



t(6-2 e -n) = -i(6 • 2 e - 1 - 2±i) - t(6 • 2 6 " 1 - " ' 



2 

e. 



- _(_i)e-i( a (!i±l) + a (!L y i)) = (-i)« a (n). 

This completes the proof. □ 

11 



3.2 A few other formulae 
Proposition 3.1 (i) We have 

s(2 e+1 + n) = s(2 e + n) + s(n) 

for < n < 2 e , (see the remark by T. Tokita concerning the Stern- sequence 
A2487 in $}). 
(ii) We have 

t(2 e+1 + n) + t(2 e + n) = (-l) e+1 s(n) 

forO<n< 2 e . 

The formulae of Proposition 13.11 have the following conjectural generali- 
sation, analogous to Conjecture 11.31 

Conjecture 3.2 (i) The series 

A{z) = E£*('(2 + n) - «(1 + n))s" =1 _ 2z + 2z 2_ 4/ + 4/ + 2^ + . . . 

l^n=0 S \ n ) zn 

satisfies 

oo oo 

^(s(2 e+1 + n) - s{2 e + n))z n = A{z T ) ^ s(n)z n 

n=0 n=0 

for all e € N. 

Similarly, the series 

J2T n(t(2 + n) + t(l + n))z n ,, , fi 7 

B(z) = _ ^»=o^ 7^ J\ =l-2z-2z 2 + 4z 3 + 6z 6 -6z 7 + ... 

sato/ies 

oo oo 

(-l) e+1 ^(t(2 e+1 + n) + t(2 e + n))z n = B{z 2 *) ^ s(n)z n 

n=0 n=0 

for all e € N. 

Proof of Proposition 13. ll The formulae hold for n = and a € {0, 1}. 

The induction step is an easy computation for odd a and obvious for even 

a. □ 

Proposition 3.3 (i) We have 

s{n) = -s{n - 2 e ) + s(n - 2 • 2 e ) + 2s(n - 3 • 2 e ), 2 e+2 < n < 2 e+3 - 2 e 

for e > 0. 

(ii) We have 

t(n) = t(n - 2 e ) - t(n - 2 e+1 ) 

for2 e+2 <n< 2 e+3 . 



12 



Proof The case e = implies n £ {4, 5, 6, 7} in assertion (i) and we have 

s(4) = 1 = -2 + 1 + 2-1 = -s(3) + s(2) + 2s(l) 
s (5) = 3 = -1 + 2 + 2-1 = -s(4) + s(3) + 2s(2) 
s(6) = 2 = -3 + 1 + 2-2 = -a(5) + s(4) + 2s(3) 
s (7) = 3 = -2 + 3 + 2-1 = -a(6) + a(5) + 2s(4) 

The induction step for e > is easy if n is even and involves the usual 

identity s(n) = s((n - l)/2) + s((n + l)/2) if n is odd. 

The proof of assertion (ii) is similar. □ 

The following result gives a few partial sums associated to the Stern 

sequence and its twist: 

Proposition 3.4 We have 

v- , x 3 e + l 



n=l 

2 C 



2 



E(-!) ns W = ^A e > X 



n=l 

2 



E*(») = ^V^^o 



2 C 



2 



EC- 1 )"*^) = 3+ i 1)e >e>0 



n=l 

Proof The first equality holds for e = and 

2 e + l 2 e 2 e 



E'W E^ + E^ 2 " -1 ) 

ra=0 n=l 

2 e 2 e 

X>(n) + E(s(n-l) + s(n)) 



n=0 n=0 n=l 

2 e 2 e 



n=0 n=l 



/ ^s x — ^ / % 3+1 3 +1 

- -#) + 3^ S (") = -l + 3— = 

n=0 

by induction. 

For the next identity one finds similarly 

oe+1 oe 

E(-i) n *w = i - E *(») = i - ^f 1 = ^ • 

n=l n=l 

The computations for the partial sums involving t(n) and (— l) n t(n) are 
analogous. □ 

We end this section with a list of a few more identities. 
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Proposition 3.5 We have 

s(3-2 e + n) = s{3-2 e - n) (3) 

for all e, n such that < e < 2", 

s{3 ■ 2 e + n) = s{3 ■ 2 e ~ l + n) + 2s{n) (4) 

for all e, n such that < n < 2 e_1 , 

t(2 e + n) = t(2 e + n- 2 e ~ 2 ) - t(2 e + n- 2 6 " 1 ) (5) 

for all e, n such that e > 2 and 1 < n < 2 e , 

s (2 e + n) = (-l) e t(2 e + n) + 2s(n) (6) 

/or a// e, n suc/i t/iat < n < 2 e+1 , 

s(2 e + n) = (-l) e t(2 e - n) - 3s(n) (7) 

/or a// e,n such that < n < 2 e ~ l , 

s (2 e - n) = (-l) e t(2 e - n) + 2s(n) (8) 

for all e,n such that < n < 2 e ~ l , 

s{2 e - n) = {-l) e t{2 e + n) + s(n) (9) 

for all e, n such that < n < 2 e . 

Proofs are easy and left to the reader. 

4 Proofs related to factorisations 

Proof of Theorem rO] We set i/j e = z(l + z 2e )(l + z + z 2 ) ITn^oi 1 ~ ^ + 
z 2 ) e_1_n . Iterating the trivial identity 

(1 + z n + z 2n ){\ -z n + z 2n ) = (1 + z 2n + z 4n ) 

we get the equivalent expression 

iP e = z(l + z 2e )l[(l + z 2n +z 2n+1 ). 
n=0 

The proof of the identity ip e = (— l) e X^n=o*(^ ' ^ e + n ) z ™ is by induction 
on e. It holds for e = 0. The induction step follows from the recursive 
definition of the sequence t(0),t(l), . . . and from the equality V'e+iX-z) = 

(i + l + z)Ve(^ 2 ). ' □ 
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Proof of Theorem II. 7t Using the Carlitz factorisation 

oo 

s( Z ) = n (i + *** + z 2n+r 

n=0 

of S(z) = Y^n=o s ( n + ^) zTl we have 



d ^2 n z 2n - 1 + 2 nA 



H(z) = -log(S(z)) = Y J 



1 2"+ 1 -l 



n=0 

The summand of index n = yields 1 !"^?^ and the sum X^Li ' ' ' can ^ e 
rewritten as 2zH(z 2 ). 

The proof of 2— regularity of i?(-z) is an easy consequence of the func- 
tional equation for H, see Theorem 11.81 below. 

Uniqueness of H defined by the functional equation H(z) = jj^r^ + 
2zH(z 2 ) follows from the fact that the map 

has a unique attracting fixpoint for formal power series (with respect to the 
obvious topology given by coefficent-wise convergency) . □ 

Proof of Theorem 11.81 We assume first that no linear form L\, . . . , L^ 
involves coefficients of degree > k. For i = 0, . . . , k — 1, we denote as before 
by p(i) the linear map 



p(i) ^2 a(n)z n =^2a(i + nk)i 



\n=0 / n=0 

Given solutions Ux, . . . , Ud, we consider a finitely generated vector space or 
module V containing Ux, . . . , U^ and the k— kernel of Ax, . . . , Ad- We have 
then 

P (i)U j =p(i)A 3 + {[x i ]L j )(U l ,...,U d ) 

where [x*]Lj € R[xx, . . . , x<j] is the linear form obtained from Lj by consid- 
ering the coefficients of z 3 . The power series p(i)Uj € V is thus a linear 
combination of Ux, . . . , U\ and of the k— kernel of Aj. The set V contains 
thus the k— kernel of Uj and Ux, ■ ■ ■ ,Ud are all k— regular. 

If there are linear forms among Lx, ■ ■ ■ ,Ld which are of degree > k, 
we introduce the k— regular series Am = zAx, . . . ,A 2 d = zAj, the series 
U d +x = zUx, ■ ■ ■ ,U 2 d = zU d , the linear forms L d+ x = zLx, ...,L 2 d = zL d G 
zi?[z][xi, . . . , Xd] and set oid+x = . . . = a 2 d = 0. We have then the equations 

U(z) = Mz) + zU{Ux{z k ), . . . , U n (z k )) 

and the identities on = U(0) for i = 1, . . . , 2d. Modifying a linear form Lj 
of degree > k by substituting all occurences of z k+l Uj(z k ) with z l Ud+j{z k ) 
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we construct an equivalent system with strictly smaller maximal degree for 
the linear forms Lj, . . . , L2d- Iteration of this construction leads eventually 
to a system containing only linear forms of degree strictly smaller than k. 

Existence and unicity of the solution follow from unicity of the attracting 
fixpoint of the dynamical system defined by the map 

Ui(z) .— > Mz) + LiCE/iC**), . . . , U d (z k )),i = 1, . . . ,d 

starting from the point («i, . . . , a<j). □ 

Proof of Theorem 11.101 The first part follows from Theorem 11.81 ap- 
plied to the identity A(z) = P{z)A{z k ). We present here however a second, 
independent proof. 

Working over the field of complex numbers and using the fact that prod- 
ucts of two /c— regular series are k— regular (cf. Theorem 16.4.1 in [2J), it is 
enough to prove the result for polynomials of degree 1. We can thus assume 
that P(z) = 1 + Xz. The coefficient of z n in A(z) = l\™=o( l + A -0 is then 
given by zero if the fc— ary expansion of n involves digits greater than 1 and 
it is given by X a otherwise, where a equals the number of ones in the k— ary 
expansion of n. This implies k— regularity of A(z). 
We have 

o/ n d , ,a, ^ ^P'(z kn )k n z k "- 1 

ra=0 v ' 

P'(z) 

The summand of index n = yields p }J and the remaining summation 

X^^Li ■ ■ ■ can De rewritten as kz k ~ 1 B(z k ). This shows that B(z) satisfies 
the functional equation 

B{z) = ^ + kz k - l B{z k ). 

p'( z ) 
Using Theorem 16.4.3 of |2J we see that the rational fraction p } J is k— regular 

if and only if all zeroes of P{z) are roots of unity (ie. if P(z) divides (z N — 1) N 
for some integer N). Theorem 11.81 implies then A;— regularity of B(z). □ 

5 Proof of Theorem 11.131 and other results involv- 
ing matrices 

Proof of Theorem 11.131 The trivial identities 



det(M(2n)) = det 
= det 



s(2n) s(2n + l) 

t(2n) t(2n + 1) 

s(n) s(n) + s(n + 1) 

-t{n) -t(n)-t(n+l) 
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and 

det(M(2ra - 1)) = det 



/ s(n) s(n + 1) 
\t(n) t(n + l) 

■det(M(n)) 

s(2n - 1) s(2n) 



i(2n - 1) t(2n) 



, ( s(n — 1) + s(n) s(n) 
V -t(n - 1) - i(n) -i(ra) 

-rW f s ( n_1 ) s ( n ) 
V t(n - 1) f (n) 

= - det(M (n - 1)) 

imply the result. □ 

Proof of Proposition 11.141 The reduction modulo 2 of the Stern se- 
quence s(0),s(l), . . . is the 3— periodic sequence 0,1,1,0,1,1,.... Indeed, 
this holds for so = 0,si = S2 = l and the recursive formulae 

s(6n) = s(3n) 

a(6n + l) = s(3n) + s(3n + l) 

s(6n + 2) = a(3n+l) 

s(6n + 3) = s(3n + l) + s(3n + 2) 

s(6n + 4) = s(3n + 2) 

s(6n + 5) = s(3n + 2) + s(3n + 3) 

imply the 3— periodicity of s(n) (mod 2) by induction. The reduction mod- 
ulo 2 of twisted Stern sequence t(0),t(l), . . . coincides with the reduction 
modulo 2 of the Stern sequence. □ 

5.1 Other results involving matrices 

The proofs of the following results are easy and omitted. 

Proposition 5.1 (i) The matrices 

s(n) s(n + 1) 

s(2 e +n) s(2 e + n + l) 

have determinant —1 for n such that < n < 2 e and determinant 1 for n 
such that 2 e <n< 2 e+l . 
(ii) The matrices 

s(n) s(n + 1) 

t(2 e + n) t(2 e + n + l) 
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have determinant (— l) e+1 forn such that < n < 2 e and determinant (— l) 6 
for n such that 2 e < n < 2 e+2 . 
(Hi) The matrices 

t(n) t(n + l) 

s(2 e +n) s(2 e + n + l) 

have determinant (— l) e+1 for n such that 2 e+1 < n < 5 • 2 e . 
(iv) The matrices 

t(n) t(n+l) 

t{2 e + n) t(2 e + n + l) 

have determinant 1 for n such that 2 e ~ 2 < n < 2 e or 7 ■ 2 e < n < 2 6+3 and 
determinant —1 for n such that 2 e < n < 7 ■ 2 e . 

6 Proof for Theorem 11.151 

For odd n we have 

s{n) = s((n - l)/2) + s((n + l)/2) = s(n - 1) + s(n + 1) . 

We have thus s( - n ~ ' s{ - n — I = \ = i -|_ 2u 2 (n) since ^2( re ) = if n is odd. 
For n even we have by induction 

s(n — 1) + s(n + 1) 



s(n) 
g((n - 2)/2) + s(n/2) + s(n/2) + g((n + 2)/2) 

s{n/2) 
s(n/2-l) + S (n/2 + l) 
s{n/2) 
= l + 2v 2 (n/2) + 2 = l + 2v 2 (n) . 

This ends the proof of assertion (i) . 

For the twisted Stern sequence we use the analogous identities 

t(n - 1) + t(n + 1) = t(n), n odd, 

t(n - 1) + t(n + 1) = -(t(n/2 - 1) + t(nj2 + 1) + 2t(n)), n even, n > 4. 

This implies assertion (ii) by checking the initial cases and the case of n S 

3 • 2 N . □ 
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